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Defects and multistability in eutectic solidification patterns 
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Abstract. - We use three-dimensional phase- field simulations to investigate the dynamics of the 
two-phase composite patterns formed upon during solidification of eutectic alloys. Besides the 
spatially periodic lamellar and rod patterns that have been widely studied, we find that there is a 
large number of additional steady- state patterns which exhibit stable defects. The defect density 
can be so high that the pattern is completely disordered, and that the distinction between lamellar 
and rod patterns is blurred. As a consequence, the transition from lamellae to rods is not sharp, 
but extends over a finite range of compositions and exhibits strong hysteresis. Our findings are in 
good agreement with experiments. 
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The spontaneous emergence of solidification microstruc- 
tures during the freezing of pure substances and alloys is 
a classic example for pattern formation outside of equilib- 
rium. These structures are of considerable practical im- 
portance in metallurgy because of their influence on the 
properties of the finished material pQ . Solidification exper- 
iments also provide a particularly well-controlled setting 
to study fundamental questions of nonlinear pattern dy- 
namics [2]. 

Eutectic alloys exhibit a triple point in their phase dia- 
gram at a temperature Te, where liquid of composition Ce 
can coexist with two solid phases a and f3 of different com- 
positions C a and C/3 . Growth from a liquid of composition 
sufficiently close to Ce leads to the formation of a compos- 
ite solid, with two basic growth morphologies: alternating 
parallel platelets (lamellae) of the two solid phases, or fi- 
bres (rods) of one phase inside a matrix of the other, with 
the centres of the fibres located on a triangular lattice. In 
cross-sections of the finished material, the two morpholo- 
gies look like stripes and round dots, respectively. This is 
reminiscent of the stripe and hexagon patterns that ap- 
pear in many two-dimensional or quasi-two-dimensional 
systems, both out of equilibrium (e. g. Turing patterns, 
convection, reaction-diffusion systems) and in equilibrium 
(e. g. thin magnetic films, block copolymers). This is of 
course not accidental: like in these systems, two antago- 
nistic effects (here, solute redistribution by diffusion and 
surface tension) lead to the emergence of a characteristic 



length scale. 

While it has been known for a long time that rods are 
preferred over lamellae when the volume fraction of one of 
the two phases is small, the precise conditions for the ap- 
pearance of these two morphologies, as well as the nature 
of the transition between them, have remained unclear. 
Here, we show by quantitative three-dimensional phase- 
field simulations of eutectic solidification that the transi- 
tion from rods to lamellae has some features that are not 
observed in other systems. We observe that regular arrays 
of rods exhibit, when their spacing is increased, a bifurca- 
tion from circular rods to symmetry-broken rods that have 
an oval or even dumbbell-like shape. The resulting stable 
arrays look like broken lamellae. There is thus a whole 
family of spatially periodic states "in between" rods and 
lamellae. Furthermore, we find that besides these periodic 
states, a large number of disordered steady- state configu- 
rations exist which consist of a mixture of rods and lamel- 
lae. Whereas such disordered states have also been found 
in other systems, they seem to be particularly stable in eu- 
tectic solidification. As a consequence of the existence of 
these numerous intermediate stable states, the transition 
from rods to lamellae is strongly hysteretic. All of these 
observations compare favourably to recent experimental 
findings. 

Let us start out by a brief review of what is known on the 
two types of periodic patterns in eutectics. We consider 
directional solidification: a sample is pulled with fixed ve- 
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locity V in an externally imposed temperature gradient G. 
The liquid freezes into two distinct solid phases of different 
compositions. If all crystallographic effects are neglected 
and the isotherms are assumed to be planar, the problem 
is invariant by rotations and translations in the plane of 
the isotherms. Under these conditions, Jackson and Hunt 
[3] obtained two families of approximate steady-state solu- 
tions for perfectly periodic lamellae and rods, and related 
the average front temperature to the lamella or rod spac- 
ing A. This front temperature is determined by the inter- 
play between the diffusive redistribution of alloy compo- 
nents through the liquid and the effect of surface tension, 
which shifts the interface temperature by an amount that 
is proportional to the local interface curvature. For the 
slow growth regime that is relevant for most experiments 
(small Peclet number, XV/D <C 1, where D is the solute 
diffusion coefficient in the liquid), the global front under- 
cooling with respect to the eutectic temperature is min- 
imal for a characteristic spacing Xjh = KV~ X I 2 , where 
the constant K depends on the thermophysical properties 
of the alloy, the volume fractions of the phases, and the 
type of pattern (lamellae or rods). 

In order to be observable in experiments, the steady- 
state solutions must be stable. Some possible instabilities 
were already discussed by Jackson and Hunt [3]. Subse- 
quent research has shown both through experiments and 
simulations that stable lamellar growth occurs within a 
range of lamellar spacings around Xjh both in thin [4H7] 
as well as in extended samples [8[[9]; for a more detailed 
review, see [TU] , Much less is known about the stability of 
rods. Recent advances have been obtained by using sam- 
ples of transparent organic alloys [TTJ[T2] . Numerically, a 
few examples of rod structures have been calculated by 
phase-field models [l3j[T4] , but to our knowledge no de- 
tailed quantitative investigation of the steady-state solu- 
tions was previously carried out. We therefore start our 
investigations by a detailed survey of rod stability. 

We use an efficient phase-field model recently intro- 
duced to simulate coupled eutectic growth (15, 16 , and 
used previously to investigate the nature of pattern insta- 
bilities of lamellar arrays in massive samples [8] as well 
as the influence of temperature profiles on the solidified 
patterns [T7] . The model has proved to be capable of pro- 
ducing results in excellent agreement with experimental 
observations [8j[T6j[T7]. The equations of the model as 
well as the method to choose appropriate computational 
parameters are discussed in detail in Refs. [8j|T6]. The es- 
sential point is that in this model, the thickness W of the 
diffuse interfaces can be chosen freely without introducing 
computational artefacts as long as it remains about an 
order of magnitude smaller that the size of the smallest 
features in the pattern. Since we are interested here in 
general aspects of eutectic solidification, we use a eutectic 
alloy with a symmetric phase diagram, which yields the 
best computational performance. In this case, all relevant 
physical parameters can be grouped to yield three charac- 
teristic length scales: the thermal length It = |m|AC/C, 
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Fig. 1: Sketch of the simulation geometry. A perfectly periodic 
rod array of horizontal spacing A and distance a between rows 
is simulated in the reduced simulation cell indicated by the 
dashed rectangle, with reflection boundary conditions. 



the diffusion length Id = D/V, and the capillary length 
d = jTE/(L\m\AC), with m the liquidus slope, AC the 
difference in composition between the two solids, 7 the 
solid-liquid surface tension, Te the eutectic temperature, 
and L the latent heat of melting. For all our simulations, 
we use Id /do = 10 3 and It /Id — 4. 

The simulation setup is similar to the one used in 
Ref. [8] . A rectangular simulation cell is filled by the solid 
composite up to the solid-liquid interface, and by the liq- 
uid phase beyond it. The isotherms (and hence also the 
solid-liquid interface) advance at constant speed V, and 
the simulation cell is moved whenever the front has ad- 
vanced by a certain amount, so that the interface is ap- 
proximately stationary in the cell. An efficient multi-grid 
scheme permits to calculate the concentration field C(x, t) 
up to a very large distance from the interface. An impor- 
tant parameter is the initial composition Co of the liquid, 
because it determines the equilibrium volume fraction 77 
of the f3 phase through the lever rule, 77 = (Co — C a )/AC. 

Since it turns out that distortions of the rod array play 
an important role for the interpretation of our results, it 
is useful to give some details about the geometry of our 
simulation cell, shown in Fig. [TJ We consider perfectly 
periodic arrays of rods, where the distance between rods 
in a "horizontal" row (along the x direction) is A, and the 
distance between rows is a. Steady states of such periodic 
rod arrays can be efficiently simulated by taking advan- 
tage of the two orthogonal planes of mirror symmetry that 
run through the center of each rod. We use the reduced 
simulation box indicated by a dashed rectangle, with re- 
flection boundary conditions on all lateral boundaries for 
all fields. The simulation cell thus contains only two quar- 
ters of rods in two corners. The picture of a full array can 
be reconstructed by adding suitably reflected and shifted 
copies of the simulation cell. The simulation is initialised 
with two circular quarter-rods, where the radius is chosen 
such that the volume fractions of the two solid phases are 
consistent with the global alloy composition. The evolu- 
tion of this initial state is then followed until a steady state 
is reached. 
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Fig. 2: Front undercooling A versus reduced spacing A = 
X/Xjh for rod arrays of two different geometries and volume 
fraction of the f3 phase of 77 = 0.3. Circles: a/X = 14/16; 
squares: a/X = 13/16. Insets: reconstructed steady states. 



A perfect triangular lattice would correspond to a/X = 
v/3/2. However, since our finite- difference discretization 
uses a regular cubic grid, our simulation box has a ratio- 
nal aspect ratio, which results in a weakly distorted array. 
In such an array, the spacings between rods in a horizontal 
row and between rods in two distinct rows are slightly dif- 
ferent. We will use the spacing A within a horizontal row 
in our subsequent plots. To investigate various spacings 
at fixed aspect ratio, we keep the size of the simulation 
box in terms of grid points as well as the ratio of the grid 
spacing Ax and interface thickness W constant (we use 
Ax/W = 0.8), whereas we vary the ratio W/d, which cor- 
responds to changing the physical size of the simulation 
box. Typically, we use 32 grid points along the smaller 
dimension of the box (of length equal to A/2, see Fig. [1]), 
which corresponds to a resolution of X/W = 51.2, suffi- 
cient to guarantee well-converged results [SJ[T6] . 

We consider first the case of a sample with volume frac- 
tion r] = 0.3 of the f3 phase (for our symmetric phase 
diagram, the eutectic point is at 17 = 0.5). In Fig. [2j we 
plot the dimensionless undercooling 

a Te — Tfront 
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as a function of the reduced spacing A = A/Aj#, where 
A j# is the minimum undercooling spacing calculated from 
the Jackson-Hunt theory for rods. Data for two differ- 
ent distorted arrays are shown. For an aspect ratio of 
a/X = 14/16, which is slightly larger than \/3/2, the un- 
dercooling curve exhibits a minimum at a spacing very 
close to the one predicted by the Jackson-Hunt calculation, 
and follows the expected behaviour A = A m i n (A + 1/A) 
[3] up to spacings of about A = 1.3. However, as the spac- 
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Fig. 3: Bifurcation diagram of the shape transition for a/X = 
13/16 and 77 = 0.3. The parameter ip is plotted versus the 
initial rod spacing. Squares and circles denote runs with dif- 
ferent initial conditions: circles, cylindrical rods; squares, rods 
elongated towards second neighbours. The top and bottom 
branches correspond to rods with a cross-section elongated to- 
wards first and second neighbours, respectively, as shown in 
the pictures of reconstructed arrays (both for A = 1.4). The 
centres of the elongated rods remain on the initial rod lattice. 

ing increases, we observe a change in shape from circular 
rods to elongated rods with the long axis oriented towards 
the second nearest neighbours (see Fig. [3] for the view of a 
reconstructed array). At the same time, the undercooling 
starts to differ markedly from the Jackson-Hunt predic- 
tion. For a slightly lower aspect ratio of a/X = 13/16, 
smaller than V3/2, the undercooling curve again exhibits 
a minimum for a slightly larger spacing, but upon increas- 
ing the spacing, this time the rods elongate in the direction 
of the first neighbours. 

In fact, both of these elongated states can be reached for 
each aspect ratio by starting from initial conditions which 
favour one direction or the other, such as elliptic rods elon- 
gated towards the first or the second nearest neighbours 
instead of simple circular rods. Thus, the shape transition 
is a bifurcation. For a more quantitative description, we 
define a shape parameter 
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with d\ and d 2 being the diameter of the rod along 
axes oriented towards the first and second nearest neigh- 
bours, respectively. A perfectly circular rod corresponds 
to i/j = 0, whilst for a deformation towards first (sec- 
ond) nearest-neighbours ip has a positive (negative) sign. 
This parameter is plotted versus the spacing in Fig. [3] for 
a/X = 13/16. Note that the bifurcation is asymmetric be- 
cause stretching of the rods towards first and second near- 
est neighbours is not equivalent. Furthermore, for this 
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Fig. 4: Comparison of the front undercooling A of rods with 
a/X = 14/16 (circles), and lamellae of half the original rod 
spacing (diamonds). In order to achieve a direct comparison 
of corresponding states, for the lamellae the reduced spacing 
was calculated as A = 2A/Aj#, where Xjh is the rod minimum 
undercooling spacing. 



aspect ratio, the single branch of weakly distorted rods 
that exists below the bifurcation threshold connects con- 
tinuously to the upper branch of solutions (horizontally 
elongated rods), whereas the second branch appears by a 
discontinuous jump of the shape parameter. The oppo- 
site behaviour is observed in the bifurcation diagram for 
the aspect ratio a/X = 14/16 (not shown): here, the sin- 
gle branch exhibits slightly negative values of ip when it 
approaches the bifurcation threshold, and connects contin- 
uously to the lower solution branch (vertically elongated 
rods). Thus, both elongated rod states are stable steady- 
state solutions, but their respective basins of attraction 
depend on the distortion of the initial rod array. Elonga- 
tion is always favoured in the direction that is stretched 
with respect to a symmetric hexagonal array. 

As can be appreciated from Fig.[3j the periodic arrays of 
elongated rods become similar, for both orientations, to a 
set of "broken" parallel lamellae, with a spacing Xi = A/2 
for the case of i/j < and Xi = a for ip > 0. Indeed, this 
is not just a geometrical similarity: in Fig. HJ it can be 
seen that the undercooling of a broken lamellar state with 
ip < is very close to the one for perfect lamellae of half 
the original spacing. Note that we have plotted the curve 
for lamellae as a function of 2 A in order to obtain a mean- 
ingful comparison in this graph. The shape change from 
circular rods to elongated rods takes place approximately 
at the intersection of the two Jackson-Hunt curves. This 
observation provides a simple criterion to predict the oc- 
currence of the morphology transition, although it gives 
no direct explanation for the mechanism that is at work. 
Note that we have not observed any reconnection of the 
"broken lamellae" up to the largest spacings investigated 
(A = 1.55). 

The behaviour of the horizontal broken lamellae in the 
case a/X = 13/16 is quite different. The undercoolings 
for lamellae of spacing a are not very different from the 
original Jackson-Hunt curve for rods of spacing A; the un- 
dercoolings of the horizontal broken lamellae thus do not 
exhibit a dramatic change with the spacing. However, 




Fig. 5: Lamellar array with a single defect that remains stable. 
The snapshot shown is a steady state. 



this branch of symmetry-broken solution terminates for A 
slightly above 1.5. The dumbbell-shaped rods shown in 
Fig. [2] split in two and reconnect to form vertical broken 
lamellae. 

In summary, the "broken lamellae" or "elongated rod" 
states have properties that are very close to the ones of 
standard lamellae. Therefore, the morphological change 
can be interpreted as a transition to a lamellar array that 
is not completed, since neighbouring rods do not connect. 
Clearly, a barrier prevents the establishment of a simple 
lamellar pattern. Further simulations reveal that the ex- 
istence of such barriers is a general property: whenever 
defects exist in a eutectic solidification pattern, they are 
very difficult to heal. A typical example is shown in Fig. [5] 
The initial condition for this simulation is an array of par- 
allel lamellae, where a "cut" has been made in the central 
lamella, thus creating two lamella terminations facing each 
other. The pattern could come back to a perfect lamellar 
array by just reconnecting the two terminations; instead, 
after an initial transient where the two endings slightly 
adapt in shape and length, the simulation evolves towards 
the steady state shown in Fig. [5j It is thus clear that there 
exists a large number of steady states that exhibit stable 
defects. Indeed, by repeating the operation of "cutting" a 
lamella, an arbitrary number of defects can be introduced 
in a lamellar pattern. Inversely, if in a perfect rod pattern 
two rods are "connected" to form a piece of lamella, the 
latter persists as a stable steady state. It is even possi- 
ble to create stable patterns which are a complete mix of 
lamellae and rods, such as shown in Fig. [6j 

The presence of multiple steady states separated by 
"barriers" implies that a system started from random ini- 
tial conditions will in general not reach a spatially periodic 
steady state. Fig. [71 shows examples of pattern evolution 
in an extended system using as initial condition a solid 
consisting of a random distribution of the two phases. For 
the eutectic composition, the emerging pattern is a com- 
plex network of lamellae, the evolution of which is very 
fast at first: reconnections between lamellae occur which 
lead towards a more ordered system. However, later on 
the evolution becomes very slow, since it is increasingly 
difficult for the system to overcome the barriers which 
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Fig. 6: Coexistence of rods and lamellae. The snapshot shown 
is a steady state. 



hinder reconnection. Interestingly, the resulting pattern 
shows, in addition to the emerging lamellar pattern, also 
a few isolated rods, which remain stable until the end of 
the simulation. For a j3 volume fraction of 77 = 0.3, mostly 
rods are formed that are disordered and elongated at first 
and then rearrange into a more ordered pattern. Note, 
however, that a few isolated "pieces of lamellae" remain. 
A series of similar simulations for various volume fractions 
shows that "pure" rod patterns are formed only for volume 
fractions below 77 « 0.25. It should be noted in passing 
that the final state for 77 = 0.3 mostly consists of rods, al- 
though it can be seen from Fig.H]that the global minimum 
undercooling is lower for lamellae at this composition. It 
is therefore clear that the criterion of lowest undercooling, 
often used in the metallurgical literature, is not sufficient 
to predict the emerging morphology. 

These findings imply that the switching between lamel- 
lar and rod patterns which is observed in real samples can 
only be explained by the presence of a sufficiently large 
driving force. In fig. [8] we explore one possibility for such 
global forcing: a change of composition with time. The 
simulation is started from a state of highly unstable paral- 
lel lamellae, which rapidly break up and form the irregular 
pattern of rods shown in the first snapshot. The average 
concentration in the liquid far ahead of the growth front 
is changed linearly with time, which implies a change in 
volume fraction from 77 = 0.2 to 77 = 0.8 in the course 
of the simulation. The figure shows that a transition be- 
tween rod and lamellar patterns can be triggered, but also 
shows that there is a strong hysteresis: the transition from 
rods to lamellae occurs at a volume fraction of 77 « 0.5, 
whilst the one from lamellae to rods occurs at 77 ~ 0.7. If 
the composition was decreased again starting from the last 
picture, the transition from rods to lamellae would happen 
again at 77 w 0.5, which shows that the patterns found for 
volume fractions between these two limits depend on the 
history of the system. 

It is interesting to notice that both the lamellar and rod 
patterns that emerge after the transitions are clearly more 
regular than the ones obtained from random initial condi- 
tions. This is a result of a non-trivial spatial organisation 
of these transitions. The formation of lamellae occurs by 




Fig. 7: (Colour online) Simulations of the symmetric model 
alloy at the eutectic composition (77 = 0.5, top row) and at 
77 = 0.3 (bottom row) in an extended system. The initial state 
is a random distribution of patches of the two solid phases. 
We have used d/W — 0.06218, and the lateral system size 
is 160 W, which corresponds to 1.6 Id and 6.5 Xjh (at the 
eutectic composition). Snapshots for both simulations were 
taken at times £/£d = 0, 0.90 and 21.87, where to — Id/D 
is the diffusion time; the last frame reflects a total solidified 
length of 127.39 A JH . 



successive connections of rods, which could lead to a ran- 
dom lamellar pattern. However, as was discussed above, if 
the initial rod array is not perfectly hexagonal, but slightly 
distorted, the shape bifurcation favours a well-defined di- 
rection of elongation, which leads to the formation of large 
patches of parallel lamellae. The transition from lamel- 
lae to rods, in turn, starts at lamella terminations and 
propagates along lamellae by successive pinch-offs of new 
rods from the ending of the lamella; this leads then to 
the breakup of neighbouring lamellae. It would be highly 
interesting to study the dynamics of these transitions in 
more detail, but this is outside the scope of the present 
work. 

All our results are in excellent agreement with experi- 
mental observations. Rods that are elongated (both to- 
wards first and second neighbours) have been observed in 
Al-Cu [18], mixed states of lamellae and rods have been 
found both in Al-Cu [19] and in transparent organic alloys 
[9]. In Fig. 7c of Ref. [20], alternating lamellae and rows 
of rods can be seen, which bear a striking similarity to our 
Fig. [6] Furthermore, it was shown that labyrinth states 
evolve indeed very slowly in time in the absence of exter- 
nal forcings [17] . Patterns of coexisting rods and lamellae 
were also observed in the peritectic Cu-Sn alloy [2T] . 

Both aspects described above - shape instabilities of 
rods and disordered and "mixed" states - have been ob- 
served in other pattern- forming systems, but there are im- 
portant differences between our findings and the ones re- 
ported in the literature. For instance, instabilities leading 
to dumbbell-shaped domains have been observed in fer- 
rofluid droplets [22], and epitaxial islands [23], and have 
been analysed theoretically in systems with a Coulomb in- 
teraction [23]. However, in all these systems a long-range 
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Snapshots are taken at t/tD = 3.2, 8.3, 10.2, 15.4, 16.7, and 
19.2; the lateral system size is 1.6 Id, the same as in Fig. [71 
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